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Abstract 



43; 

Interpretation of long-range rapidity correlations in terms of the fluctuating 
|^-i rapidity density distribution of the system created in high-energy collisions 

is proposed. When applied to recent data of the STAR coll., it shows a 
substantial asymmetric component in the shape of this system in central 
Au-Au collisions, implying that boost invariance is violated on the event-by- 
event basis even at central rapidity. This effect may seriously influence the 
hydrodynamic expansion of the system. 
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1. It is now widely recognized that long-range correlations (LRC) in ra- 
pidity originate at the early stages of the collision, before the longitudinal 
^ ■ expansion separates the particles by large distances. Such correlations can 

thus be used as a probe of the initial conditions of the evolution. This is 
particularly interesting for hydrodynamic description of particle production, 
as the initial conditions strongly influence the evolution of the system (called 
henceforth " a fireball" ) expanding according to the rules of hydrodynamics 
[l| . The event-by-event fluctuations of the initial conditions in the transverse 
plane were already shown to induce several interesting features in the trans- 
verse momentum correlations observed in the final state [2J . For example the 
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recently observed so-called triangular flow is entirely due to the asymmetric 
fluctuations in the transverse plane |3J]. In the present paper we discuss the 
effects of the initial state fluctuations on the correlations between the longitu- 
dinal momenta. We show how the measurements of the LRC in rapidity can 
be interpreted in terms of the event-by-event fluctuating, shape of the fireball. 
We argue, using the relevant data of STAR coll. |4|, |5(, that such analysis 
can uncover some hitherto unobserved features of particle production. 

A special case of LRC are forward-backward correlations where one com- 
pares particle distributions in two intervals located symmetrically in the for- 
ward and backward hemispheres. They were extensively studied since the 
early times of high-energy physics |6J. In most of these studies only the 
global density fluctuations were considered [7|, |8|, M, llO| and data were inter- 
preted as evidence for strong event-by-event fluctuations of the multiplicity 
of the produced particles. With the increasing precision of data and larger 
observed particle densities, however, it seems useful to consider the more 
general scenario, with the event-by-event fluctuations of both multiplicity 
and shape of the created system. In the present paper we show that this ap- 
proach allows to obtain some direct information about the object produced 
in a high-energy collision. 



To this end we apply the recently proposed method [ll|, |12|, |13| of a 



systematic study of the factorial moments of multiplicity distribution in sev 



eral well separated rapidity bins (see also [l4|). Let us add that, as shown 



in 



111 . |13I |. such measurements allow also to discriminate between various 
models of the multiparticle production and thus to understand better the 
mechanism of such processes. 

In the next section the problem is formulated in terms of the generating 
functions. In Section 3 the data of STAR collaboration [4| are analyzed and 
it is shown that they imply existence of a substantial asymmetric component 
in the fireball shape for central Au-Au collisions. A discussion of this result 
is given in section 4. Our conclusions are listed in the last section. 

2. Consider a fireball created in a single collision and a rapidity bin A*. 
The number of particles in Aj is a random number, whose distribution de- 
pends on the initial conditions of the collision, on parameters of the fireball 
evolution and on details of hadronization. All these factors may be summa- 
rized in a set of parameters, Q = (gi, q2, ■■■) which are also random numbers. 
They may be, for instance, the impact parameter of the collision, the number 
of participating nucleons, inhomogeneities of the expanding fluid, etc. Let us 
denote the average number of particles falling into Aj at a given Q by n^Q). 



It is related to the corresponding particle density by 

ni{Q) = / p(y)dy « p(Vi)\ (1) 

JAi 

We call dynamic the fluctuations of n« resulting from the randomness of 
(qi,q 2 , ■■■)■ Our purpose is to estimate these fluctuations. 

Even for a fixed Q, however, the actual number of particles in Aj fluc- 
tuates around rij. This is the noise we want to correct for. To this end we 
assume that it is possible to choose the parameters (<7i,<72, •••) so that for 
each set Q the fluctuations around n^Q) are dominantly random i.e. ap- 
proximately Poissonian. Note, that we do not have to know the set Q. It is 
enough to assume that it exists. Under this assumption, for B bins at given 
Q the probability distribution for the occupation numbers rii, . . . , n B is 

- n n 
P(n 1 ,...,n B ;n 1 ,...,n B )=p(n 1 ;ni)...p(n B ;n B ); p(n;n) = e~ n — (2) 

n! 

where the argument Q has been omitted. The observed distribution is the 
average over Q, or equivalently over the averages fii(Q)j. 

P(ni,...,n B )= dn 1 ..Jn B W(n 1} ...,n B )p(n 1 -,n 1 )...p(n B -,n B ) (3) 

where W{n\, ...,n B ) is the probability distribution of the set [n 1; ...n B ], char- 
acterizing the distribution of the densities of the produced fireballs, i.e. the 
basic quantity of interest. 

From the well-known property of the Poisson distribution 
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-p(n; n) = n k (4) 



(n-k)\ J 
one easily derives 

Fii,..,i B = I dn 1 ...dn B W(n 1 , ...,n B )f% , ...,n*£ = {n % i,...,fi*g) w , (5) 



lr To illustrate the idea, consider a Monte Carlo simulation in which the probability den- 
sity of finding a number of particles at some momentum is not fixed but depends on some 
(random) parameters Q. The probability density at a given Q is our p(y). Fluctuations 
resulting from fluctuations of Q are the dynamical fluctuations. The remaining random 
event-by-event fluctuations are our purely statistical fluctuations. 



where F iu __ t i B are the factorial moments of the distribution (J3J): 

\(rii-zi)! {n B -i B )\/ 

Eq.([5]) shows that measurement of factorial moments of the observed mul- 
tiplicity distribution gives directly the moments of the fluctuating fireball 
densities |l5| . 

3. As an application of ([5]) we shall show that the published results of the 
STAR collaboration [4J give evidence for a substantial asymmetric component 
in fluctuations of the fireball created in Au-Au collisions at y^i = 200 GeV. In 
this experiment the fluctuations of multiplicity observed in two rapidity bins, 
symmetric with respect to y c . m . = [at 0.8 < \y\ < 1.0], were measured for 
various centralities, selected according to the number of particles observed in 
the central bin, located also symmetrically around y c .m. = 0. At the highest 
centrality and at the highest distance between the bins, the measurements 
I, B give 

Dj f =<n)> - <n f > 2 = 350 ± 17; 
D) b =< n f n b > - <n f > 2 = 202 ± 17; < n f >= 96 ± 5. (7) 

where the indices / and b refer to the forward and backward bins, respectively, 
and rif and n& denote the actually observed numbers of particles in these bins. 
The factorial moments are related to ([7]) by 



F 20 =< n f (n f - 1) >=< n 2 f > w = D 2 ff + < n f > 2 - < n f > 



Fn =< rifn b >=< nfn b >w= D 2 fb + < nf > 2 (8) 



Noting that (n/ ± nf,) 2 = n\ + n 2 . ± 2nfn\, we thus obtain for the asymmetric 
and symmetric fluctuations 

D 2 _ = i < (n f - n b ) 2 > w = i[F 20 - F n ] = 26 ± 12; 

> 2 = ^20 + F n ]- < n f > 2 = 228 ± 12. (9) 



D + = 15.1 ±0.4. (10) 
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(n f + n b ) 2 


> w — < rif 
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we have 








D_ 


. = 5.1 ±1.2 



One sees that, although the symmetric fluctuations dominate, there is also a 
substantial asymmetric component. Indeed, the ratio D_/D + rj 1/3. Thus 
one has to conclude that created fireballs are not necessarily symmetrid_|. 
This observation implies that the standard assumption of boost invariance is 
violated on the event-by-event level even at y ~ 0. As the effect is expected to 
be stronger at the early times (because the expansion has a natural tendency 
to smooth out the original inhomogeneities), this observation may have im- 
portant consequences for the theoretical description of the process (e.g. for 
the hydro calculations). 

It is interesting to see how this result compares with pp collisions. For 
this case the STAR collaboration J4| gives Dj f = 0.572 ± 0.030 and Dj b = 
0.027 ± 0.003. Unfortunately, the value of < rif > was not given in [4|. To 



obtain a rough estimate, we have used < rif >= 0.46 ± 0.03 taken from [16 
This gives £>_ = 0.21 ± 0.05 and D + = 0.26 ± 0.04, indicating that in this 
case (i) the relative fluctuations are strongeio and (ii) the asymmetry of the 
fireballs is even more important. This was to be expected from the earlier 
observation |8[ that the UA5 pp data (l7[ are consistent with the presence 
of two asymmetric contributions (most likely representing remnants of the 
forward and backward moving projectiles). 

4. The results discussed in the previous section demand a technical com- 
ment. The point is that the STAR measurements were performed at a fixed 
number N c h of particles observed in the central bin and then averaged over 
N c h. The importance of this condition was first noted by Lappi and McLerran 
who analyzed this effect in detail [14]. Here we add only a remark that al- 
though this procedure may indeed change significantly D + , it does not affect 
D-. This can be seen from the identity 

D - = \[ D2 ff-< n f>- D lf\ (H) 

showing that D 2 _ is insensitive to the order of averaging over N t 
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The observed asymmetries find a natural explanation if, at RHIC energies, 
the remnants of the projectiles are still present even in the central rapidity 



18, 191. Indeed, the contributions from the forward and backward 



region 

moving projectiles are naturally asymmetric 18[. Since they are expected 
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Obviously, for a symmetric fireball Z?_ = 0. 

pp 



3 As is seen from the inequality [D / < n >] pp > [D/ < n >}auAv 



to fluctuate quasi-independently, they produce -generally- an asymmetric 
fireball. 

Since the STAR data give only moments of the second order, it is clear 
that introducing such two components gives enough freedom to obtain a cor- 
rect value for D-. It is not clear, however, if these two contributions are 



really sufficient to describe the physics 19|, |20j of the observed asymmetry. 
It may well be that it is necessary to add a third (approximately symmetric) 
component, describing the multiple parton-parton interactions [21]. Mea- 
surement of moments of higher order is necessary to answer this interesting 
question. The higher moments are also needed if one wants to pin down more 
precisely the details of the distributions. 

In this context one may add an obvious remark that the relevant mea- 
surements at LHC energies shall be also very interesting and actually may 
help to disentangle this problem. Indeed, at these much higher energies, the 
remnants from the projectiles in the central rapidity region are expected to 
be small and perhaps even completely die out. Therefore one expects their 
contribution to asymmetry to be significantly smaller. 

5. Some comments are in order. 

(i) As already mentioned in Section 2, the density p(y), whose fluctua- 
tions we propose to study, summarizes effects of all processes leading to the 
observed final state. In particular, it can be modified by hadronization: the 
density of gluons is replaced by densities of particles and resonances. Our 
poissonian Ansatz ([2]), ([3]) removes the purely statistical fluctuations, but 
does not remove the possible dynamical fluctuations which may be induced 
during hadronization, for example the resonance production. These fluctu- 
ations are, for obvious reasons, much more difficult to control, as they are 
model-dependent. We would like to emphasize, however, that these are the 
genuine dynamical effects which provide relevant physical information about 
the process of particle production. Therefore they should be included in the 
analysis. 

(ii) An estimate of the hadronization effects is of course of great interest, 
since otherwise it is not possible to make definite statements about the initial 
conditions which are of main interest. We have therefore worked out how the 
resonance production can modify our main results presented in Section 3. We 
have considered an extreme situation in which all observed pions are decay 
products of resonances. For p production, the correction to D- does not 
exceed 8-10%, for the transverse momentum of the p up to 1 GeV. Similar 
results are obtained for u production and for production of heavy (1.5 GeV) 

6 



clusters decaying into 3 particles. The corrections can reach up to 20% in pp 
collisions, still inside the quoted errors. The reason for this small effect is the 
small bin size Arj = 0.2, which makes it unlikely to have more than one decay 
particle in the bin. We therefore conclude that neither resonance production 
nor clustering effects can explain the asymmetry observed in data. 

(iii) It is of course not excluded that another, hitherto unknown, effect 
can be responsible for the observed asymmetry. For the moment, however, 
its most likely explanation is the asymmetry in the initial conditions. 

6. In conclusion, it is argued that the systematic study of the factorial 
moments of the multiplicity distribution in several rapidity intervals repre- 
sents a powerful tool allowing to investigate, on event-by-event basis, the 
longitudinal structure of systems ("fireballs") created in high-energy colli- 
sions. It was shown that, when applied to the data of STAR collaboration 
[4j, this method allows to uncover the importance of asymmetric fireballs 
produced in the symmetric Au-Au collisions. This result seems interesting, 
since it suggests that the hypothesis of boost invariance is violated on the 
event-by-event level even in the central rapidity region. It is also interesting 
to note that the effect seems more significant for pp collisions, thus indicating 
a violation of boost-invariance in "elementary" collisions. We feel that these 
observations should be seriously taken into account in modeling the particle 
production processes. 

It should be emphasized that the method we discuss is very general and 
flexible. It can be applied to any specific sample of events, e.g., those asso- 
ciated with a large transverse momentum jet and/or selected according to 
overall multiplicity, transverse momentum and many others. Hopefully, the 
coming measurements at LHC will be able to exploit its full capacity. 
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